Biological Conservation 57 U 991) 73-88 


- N95- 24096 


I s 


Measuring Forest Landscape Patterns in the Cascade 
Range of Oregon, USA 



William J. Ripple 

Environmental Remote Sensing Applications Laboratory, Department of Forest Resources. 
Oregon State University, Corvallis, Oregon 97331, USA 

G. A. Bradshaw 

Environmental Remote Sensing Applications Laboratory and Department of Forest Science, 
Oregon State University, Corvallis, Oregon 97331, USA 

& 



Thomas A. Spies 

USDA Forest Service, Pacific Northwest Research Station and 
Oregon State University, Corvallis, Oregon 97331. USA 

(Received 24 April 1990; revised version received and accepted 8 October 1990) 


ABSTRACT 

This paper describes the use of a set of spatial statistics to quantify the 
landscape pattern caused by the patchwork ofclearcuts made over a 1 5 -year 
period in the western Cascades of Oregon. Fifteen areas were selected at 
random to represent a diversity of landscape fragmentation patterns. 
Managed forest stands ( patches) were digitized and analysed to produce both 
tabular and mapped information describing patch size, shape, abundance and 
spacing, and matrix characteristics of a given area. In addition, a GIS 
fragmentation index was developed which was found to be sensitive to patch 
abundance and to the spatial distribution of patches. Use of the GIS-derived 
index provides an automated method of determining the level of forest 
fragmentation andean be used to facilitate spatial analysis of the landscape 
for later coordination with field and remotely sensed data. A comparison of 
the spatial statistics calculated for the two years indicates an increase in forest 
fragmentation as characterized by an increase in mean patch abundance and a 
decrease in interpatch distance, amount of interior natural forest habitat, and 

73 

Biol Conserv. 0006-3207/91/S03-50© 1991 Elsevier Science Publishers Ltd, England Pnnted 
in Great Britain 


PRW,W " M0 warn wr nun* 



74 


William J. Ripple , 0 A. Bradshaw, Thomas A. Spies 


the GIS fragmentation index. Such statistics capable of quantifying patch 
shape and spatial distribution may prove important in the evaluation of the 
changing character of interior and edge habitats for wildlife. 


INTRODUCTION 

The Douglas-fir Pseudotsuga menziesii (Mirb.) Franco forests of western 
Oregon have been cut extensively during the past 40 years. In these forests, 
clearcuts, new plantations, and second-growth stands now exist on the 
landscape formerly dominated by extensive old-growth forests and younger 
forests resulting from fire disturbance (Spies & Franklin, 1988). Con- 
sequently, the landscape has become more spatially heterogeneous. Some ol 
the effects of this newly created landscape on the forest ecosystem are 
immediately apparent. For example, the amount of old-growth forest 
habitat for interior species such as the northern spotted owl Strix 
occidentals caurina is altered with forest harvesting practices. Less well 
understood are the long-term and more subtle interactive effects on 
ecosystem processes (e.g. changes in species diversity and abundance, 
nutrient cycling and primary forest productivity), which may occur as a 
result of the changed forest mosaic. It is generally accepted that wildlife 
ecology and behavior may be strongly dependent on the nature and pattern 
of landscape elements (Forman & Godron, 1986), but few precise 
measurements relating these changes to landscape spatial alteration have 
been made. Properties of forested landscapes such as patch size, the amount 
of edge, the distance between habitat areas, and the connectedness of habitat 
patches have a direct influence on the flora and fauna (Thomas, 1979; Hams, 
1984; Franklin & Forman, 1987; Ripple & Luther, 1987). For these reasons, 
models and monitoring schemes are urgently needed for prescribing t e 
location, size, and shape of future harvest units and old-growth habitat 
patches. With proper design, these forest landscapes should be able to 
achieve desired habitat values and maintain biological diversity (Noss, 1983, 

Harris, 1984; Franklin & Forman, 1987). 

As a first step to studying forest landscape pattern, the spatial character of 
the landscape must be quantified to relate ecological processes to landscape 
configuration. Numerous methods and indices have been proposed for these 
purposes (e.g. Forman & Godron, 1986; Milne, 1988; O’Neill ai, 1988). 
Landscape pattern can be quantified using statistics in terms of the 
landscape unit itself (e.g. patch size, shape, abundance, and spacing) as well 
as the spatial relationship of the patches and matrix comprising the 
landscape (e.g. nearest-neighbor distance and amount of contiguous matrix)^ 
A selection of these measures can therefore describe the several aspects of 
fragmentation which occur as the result of forest harvesting practices. 
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Because spatial analyses are often cumbersome due to large amounts of 
data, it is desirable that the analysis be automated. Furthermore, it should be 
amenable for study at various scales in concert with a variety of data types. 
Automated systems such as geographical information systems (GIS) address 
such needs (Ripple, 1987, 1989). With the development of GIS, the ability to 
readily measure spatial characteristics of landscapes in conjunction with 
field and remotely sensed data has become possible. 

The overall goal of the present study was to test the feasibility of 
measuring forest landscape patterns accurately using a set of spatial 
statistics and a GIS. Specifically, we set out to apply a set of statistics to a 
series of forested landscapes in the Cascade Range of western Oregon to ( lj 
assess the sensitivity of these statistics to characterize landscape pattern; (2) 
develop and test a GIS-derived index to measure fragmentation; and (3) 
quantify the change and types of forest fragmentation through time in our 
study area. The results of this pilot study are intended to aid the forester, 
wildlife biologist, and land manager in assessing change in wildlife habitat 
and alteration in ecosystem processes as a result of forest fragmentation. 


STUDY AREA 

The study area consisted of the Blue River and Sweet Home ranger districts 
of the Willamette National Forest. According to Franklin and Dyrness 
(1973), these ranger districts lie primarily in the western hemlock Tsuga 
heterophylla and pacific silver fir Abies amabilis vegetation zone, with the 
major forest tree species consisting of Douglas-fir, western hemlock, pacific 
silver fir, noble fir Abies procera, and western redcedar Thuja plicaia. The 
climate can be characterized as maritime with wet, mild winters and dry, 
warm summers. Under natural conditions, Douglas-fir is the serai dominant 
at elevations below 1000 m, where it typically develops nearly pure, even- 
aged stands after fire. Large areas are covered by old-growth Douglas-fir/ 
western hemlock forests in which Douglas-firs are over 400 years old. Fires 
during the past 200 years have created a complex mosaic of relatively even- 
aged natural stands throughout the study area. Superimposed on this 
natural mosaic is a second component of pattern complexity resulting from 
timber harvesting over the last 40 years. 


METHODS 


Data acquisition 

Landscapes were classified in terms of managed and natural elements. 
Managed forest stands (typically young forest plantations of up to 
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approximately 40 years old established after clearcutting) were defined as 
internally homogeneous units (or managed patches) embedded in a matrix 
consisting of natural, uncut forest. This two-phase mosaic is a simplified 
representation of a more complex system in which many patch and matrix 
types exist. For the present analysis, however, we chose to classify the 
elements of the system as either ‘managed’ or ‘natural’. A managed patch 
consisted of one, or (if adjacent) more than one, unit that was clearcut in the 
past. 

Maps of forest patterns were constructed for the years 1972 and 1987. 
Managed forest patches for 1972 were mapped using high-altitude infrared 
images (scale 1 : 60 000) and were transferred to US Forest Service vegetation 
maps using a zoom transfer scope. The 1987 data were acquired from 
vegetation maps produced in 1987 by the Forest Service at a scale of 1 : 1 5 840 
for the entire Willamette National Forest. Fifteen forested landscapes were 
chosen at random from the Blue River and Sweet Home ranger districts' 
vegetation maps for this study. Each sample landscape consisted of a 
rectangle representing approximately 3-5 x 5 0 km (1750 ha). To qualify for 
selection, each landscape was required to be forested, lie within US Forest 
Service land, and contain at least one managed patch by 1972. Other 
landscapes were rejected during the random selection process. Several of the 
selected landscapes included roadless areas with much uncut forest land. 
The locations, sizes, and shapes within the managed forest patches and the 
natural forest matrix were digitally recorded as polygons for both the 1972 
and 1987 data. The total number of digitized patches for 1972 and 1987 were 
150 and 298, respectively. 

Landscape statistics 

The degree of fragmentation sustained by the forest matrix which 
characterizes a given landscape may be described as a function of the 
varying size, shape, spatial distribution, and density of clearcut patches 
(Burgess & Sharpe, 1981). Thus the degree of fragmentation can be 
measured in a number of ways. Because a single statistic is usually deemed 
insufficient to capture the entire spatial character of the landscape, a suite of 
statistics was selected. Five groups of statistics were used to quantify 
landscape heterogeneity and pattern for each of the fifteen areas for each of 
the two years: (1) patch size; (2) patch abundance; (3) patch shape; (4) patch 
spacing; and (5) matrix characteristics. Patch sizes and shapes were only 
determined for interior patches (i.e. patches that were not truncated by the 
borders of the landscape study sites). 

Patch size for each year and sample area were expressed in terms of the 
average patch area and average patch perimeter. The area and perimeter of 
each patch were computed using the digitizing routine mentioned above. 
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The means of the patch areas and perimeters for each landscape for each 
year were also calculated. The second set of statistics, ‘patch abundance’, 
includes a measure of the patch density (expressed as the number of 
managed forest patches present per landscape study area) and percent in 
patches (the percent of the total landscape area occupied by managed 
patches). Means of these two statistics were calculated for each landscape 
and year. Because patch shape and patch spacing statistics can involve more 
than one variable in their calculation, some background discussion of the 
equations is included. 

Patch shape was measured in three different ways: (1) the simple ratio of 
patch perimeter to patch area; (2) the fractal dimension; and (3) a diversity 
index. All three indices are a function of the perimeter and area of a given 
patch. The application of these similar measures on the same data set 
afforded an opportunity to compare their ability to detect spatial pattern. 

The fractal dimension, D, was used to quantify the complexity of the shape 
of a patch using a perimeter-area relation. Specifically 

P*A di2 (1) 

where P is the digitized patch perimeter and A is the patch area (Burrough, 
1986). The fractal dimension for each sample area and year was estimated by 
regressing the logarithm of patch area on its corresponding log-transformed 
perimeter. The appeal of fractal analysis is that it can be applied to spatial 
features over a wide variety of scales. A fractal dimension greater than 1 
indicates a departure from a euclidean geometry, i.e. an increase in shape 
complexity. As D approaches 2, the patch perimeter becomes ‘infilling’ 
(Krummel et at . , 1987). 

A similar index is the diversity index, DI, which was used to express patch 
shape as 



where the variables are defined as above in eqn (1) (Patton, 1975). 
Theoretically, the diversity index increases to 1 as the unit shape approaches 
a circle, similar to the case of the fractal dimension. However, in contrast, the 
diversity index increases without limit as patch shape becomes more 
complex. 

Patch spacing was characterized by measures of the mean nearest- 
neighbor distance and a measure of dispersion. The mean nearest-neighbor 
distance was calculated manually with a scale by measuring the distance 
from the centroid of each patch to the centroid of its nearest neighbor and 
computing the mean distance for the sample landscape. The centroid of each 
patch was determined through an ocular estimate procedure. Clark and 




Proximity to Managed Patches 
50m Cells 



Fig. 1. An example of a proximity map using the spread function in the pMAP geographic 
information system. Cell values in the matrix were assigned based on their distance to the 
nearest managed forest patch (areas shown in white). The GIS fragmentation index (GISfrag) 
was determined by calculating the mean of all grid cell values on the proximity map. 
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Evans (1954) developed a measure of dispersion of patches using the equation 

R = 2p 1;2 r (3) 

where R is dispersion, r is the mean nearest-neighbor distance, and p is the 
mean patch density (number of patches per unit area; see also Pielou, 1977, 
p. 155). Dispersion is a measure of the non-randomness of the patch 
arrangement. In a random population, /?=1; R less than 1 indicates 
aggregation of the patches, while R greater than 1 indicates that the patch 
population forms a regular dispersed pattern or spacing. 

A fifth group of statistics was calculated which predominantly reflects the 
character of the matrix (in this case, natural forested land) as opposed to 
managed patch configuration (clearcut areas): namely, a GIS-derived index 
(christened GISfrag), matrix contiguity, interior habitat, and total patch 
edge. To determine contiguity, a mylar sheet with an 8 x 8 grid consisting of 
64 cells each with a size of 27 ha was overlaid on the sample landscape maps. 
The largest contiguous natural forest area (i.e. the largest number of 
contiguous grid cells) in the sample landscape was recorded. Using this 8x8 
grid, the contiguity index could potentially range from 0 (a landscape with 
no natural forest patches greater than 27 ha, i.e. highly fragmented) to 64 (a 
landscape with no managed stands, and hence no fragmentation). A 27-ha 
cell size was chosen because it was considered to be a viable habitat patch 
size and fits within the structure of the existing landscape. 

The GIS fragmentation index (GISfrag) was computed by first producing 
a proximity map using the SPREAD function in the pMAP GIS software 
(Spatial Information Systems, 1986). This procedure assigned cell values 
based on the distance to the managed forest patches (i.e. a distance of one cell 
away from a managed patch was equal to 50 m, a distance of two cells was 
equal to 100 m, and so forth; see Fig. 1 for an example of a digital map 
illustrating all matrix distances to managed patches). A GISfrag was 
computed as the mean value of all the grid cell values on the proximity map, 
including the managed patches which were assigned values of zero. Large 
mean values reflected a low degree of forest fragmentation while maximum 
fragmentation occurred when the mean values approach zero. 

The spread function in pMAP was also used to calculate the amount of 
interior forest habitat. Interior forest was defined as the amount of natural 
forest remaining after removing an edge zone of 100 m (approximately two 
tree heights) into the natural forest matrix (Fig. 2). The mean total edge was 
simply the average of the total managed patch edges for each landscape. A 
non-parametric test (Wilcoxon rank sum) for the difference between mean 
values was performed for each landscape and year for each of the variables 
listed above to assess the ability of these variables to reflect landscape 
change. 
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Interior Forest Habitat 
50m Cells 


Symbol 


Label 


Percent 


Interior Habitat 57.73 


Edge Effect (100m) 21.87 

Managed Forest 20.40 


Fig. 2. An example of an intenor forest habitat map generated by pMAP. The interior 
forest habitat is defined here as the amount of natural forest remaining after removing an 
edge zone of 100 m into the natural forest. 
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RESULTS AND DISCUSSION 


Landscape structure 

The degree of forest fragmentation as measured by the statistics discussed 
above increased significantly over the 15 years spanning 1972 to 1987 (Table 
1 and Fig. 3). Overall, patch area and patch perimeter decreased by 17% and 


1972 



Mean Patch Perimeter (m) 

1793 

Mean Patch Area (ha) 

14.7 

Perimeter/ Area 

0.014 

Fractal 

1.26 

Diversity Index 

136 

Percent in Patches (%) 

8.5 

Patch Density (no./area) 

12 


Contiguity 

34 

Interpatch Distance (m) 

841 

GIS Fragmentation (km) 

0.5 

Dispersion 

1.40 

Interior Habitat (%) 

79.7 

Total Patch Edge (km) 

20.2 


0 km 5 

SCALE 


1987 



Mean Patch Perimeter (m) 

1387 

Mean Patch Area (ha) 

8.1 

Perimeter /Area 

0.020 

Fractal 

1.28 

Diversity Index 

1.42 

Percent in Patches (%) 

27.3 

Patch Density (no V area) 

38 


Contiguity 2 

Interpatch Distance (m) 498 

CIS Fragmentation (km) 0.2 

Dispersion 1.47 

Interior Habitat (%) 36.8 

Total Patch Edge (km) 68.0 


Fig. 3. A graphical and statistical illustration of landscape change between the years 1972 
and 1987 for one of the landscapes used in the study. 






TABLE 1 

Descriptive Statistics for 1972 and 1987 Landscapes 


Variable 

Mean 

Standard 

deviation 

M inimum 

Maximum 

p value 

Managed patch size 

Mean patch perimeter (ml 

2 075 

512 

1 599 

3492 


1972 

Mean patch perimeter (m) 

1 848 

452 

1211 

2 797 

0096 

1987 

Mean patch area (ha) 1972 

19 5 

9*1 

1-7 

42-7 


Mean patch area (ha) 1987 

16 2 

8-3 

6-6 

38-5 

0027 

Managed patch shape 

Perimeter/area 1972 

0013 

0-0037 

0007 

0-022 


Perimeter/area 1987 

0 018 

0-0053 

0-007 

0-028 

0-027 

Fractal 1972 

1 26 

0-02 

1-24 

1*28 


Fractal 1987 

1-27 

0-02 

1-24 

1-32 

0-023 

Diversity index 1972 

1-38 

0-11 

1 22 

1-53 


Diversity index 1987 

1-40 

0-09 

1-23 

1-52 

0-027 

Managed patch abundance 

Percent in patches 1972 

9-5 

7-0 

1-3 

27*3 


Percent in patches 1987 

18 2 

9-6 

1-3 

38-1 

0-001 

Patch density (no./area) 

10-3 

8-2 

2 

28 


1972 

Patch density (no./area) 

19-6 

9-7 

2 

38 

0-001 

1987 

Managed patch spacing 

Nearest-neighbor distance (m) 

928 

251 

589 

1375 


1972 

Nearest-neighbor distance (m) 

661 

247 

454 

1375 

0*001 

1987 

Dispersion 1972 

1 22 

0-27 

0-79 

1-66 

1-000 

Dispersion 1987 

1-24 

0-23 

0-79 

1-60 

Matrix characteristics 

GIS fragmentation (km) 

0-8 

0-5 

01 

1-8 


1972 

GIS fragmentation (km) 

0-5 

0-5 

01 

1-8 

0-003 

1987 

Matrix contiguity 1972 

36-6 

19 9 

1 

59 


Matrix contiguity 1987 

20-6 

19 1 

1 

59 

0-003 

Interior habitat (%) 1972 

78 2 

16-1 

38-4 

96-9 


Intenor habitat (%) 1987 

60 2 

19-7 

25-6 

96-9 

0-001 

Total patch edge (km) 1972 

19 9 

15-6 

2-8 

59 5 

0-001 

Total patch edge (km) 1987 

38-4 

19-8 

2-8 

716 


These statistics were calculated using patch and matrix variables of the fifteen landscapes to 
investigate the change in forest fragmentation over lime. A non-parametric means test 
(Wilcoxon rank sum) was performed to evaluate the statistical significance of the differences 
in the means (see p values). 
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11% or from 19-5 to 16-2 ha and 2075 to 1848 m, respectively. Although some 
1972 patches ‘grew’ by 1987 as a result of the coalescence of two or more 
managed patches, individual managed patch size on the whole decreased. 

The fractal dimension, the perimeter-to-area ratio, and diversity index 
indicate a statistically significant increase in patch complexity over time. A 
contributing factor to this increased complexity may be attributed to the 
linking of several smaller managed patches which often produces an 
irregularly shaped patch. The fractal dimension and diversity index appear 
to be fairly robust measures of an 'average' patch shape and may be capable 
of discerning subtle changes in patch configuration which are difficult to 
assess by visual inspection alone. In contrast, because it is not scale- 
invariant, the perimeter-to-area ratio must be interpreted carefully. 

Patch density increased over the 15 years by 98%. The percent of the 
landscape in managed patches nearly doubled in the 15-year period from 
9-5% to 18-2%. The mean of the index of contiguity decreased by 44% by 
1987, reflecting the increased number of clearcuts made in the area. Not 
surprisingly, the mean nearest-neighbor distance also decreased significantly 
from 928 to 661m. Dispersion, however, did not register the change in 
spatial distribution of patches over time. There was no difference in 
dispersion from 1972 to 1987 (1-22 versus 1-24), which represents a regular, 
dispersed spacing pattern for both dates. Since dispersion is a function of 





0 km 5 

SCALE 



0.7 km (14 cells) 


0-3 km (6 cells) 


0.1 km (2 cells) 


Fig. 4. The GIS fragmentation index (GISfrag) for a progression of landscape frag- 
mentation levels. Low levels of GISfrag reflect high amounts of forest fragmentation in 
terms of number of clearcut patches and amount of patch aggregation. This figure consists of 
six different landscapes and does not represent the same landscape over time. 
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both patch density and interpatch distance, we expected it to change because 
these other variables were statistically different over time. However, the 
mean change in these two variables was in opposite directions, i.e. patch 
density increased while nearest-neighbor distance decreased. We conclude 
that their combined effects canceled one another. The amount of interior 
forest habitat decreased from 78-2% in 1972 to 60-2% in 1987. This decrease 
in interior habitat (18% loss) was approximately twice the areal increase in 
managed patches (8-7%). 

The GISfrag mean decreased by 0-3 km in 1987, indicating an increase in 
forest fragmentation (Table 1). Figure 4 shows the GISfrag for a progression 
of landscape fragmentation levels. The GISfrag seems to be sensitive to the 
abundance of patches and the amount of unfragmented contiguous natural 
forest. Figure 5 demonstrates the influence of patch spacing and the amount 
of contiguous natural forest on the GISfrag by showing the GISfrag for a 
pair of landscapes that had the same patch density but a different level of 
patch aggregation. It should be noted that the GISfrag is only comparable 
among study areas that are the same size and shape, because the truncation 
of the matrix as it intersects the study area boundary may influence the value 
of this index. 


LANDSCAPE A 



Patch Density 22 

Percent in Patches 13J 

CIS Fragmentation (km) OJ 


0 km 5 

1 ■■■■■■ l 
SCALE 


LANDSCAPE B 



Patch Density 22 

Percent in Pitches 14.8 

G1S Fragmentation (km) 0.5 



Fig. 5. An example of two landscapes (A and B) having the same patch density (22) but 
differing level of patch aggregation reflected by correspondingly different GIS fragmentation 
indices. The accompanying figure to A and B is a graph of percent of each landscape versus 
distance of each point in the landscape to a clearcut (managed) forest patch. The B profile 
shows more of the landscape with greater distances to the nearest managed patch than A. 
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Ecological and management implications 

An examination of the results reveals the varying ability of each analysis to 
describe change in patch and matrix characteristics. The nature and amount 
of change detected by these different landscape statistics have significance to 
the ecology and management of forest landscapes. Changes that were 
evident in patch-level characteristics indicate a trend toward smaller units 
and a slight trend toward more irregular units. These changes may reflect the 
increased cutting on steeper, more irregular terrain, more careful ‘fitting’ of 
the cutting units as the available cutting area decreased, an effort to optimize 
big game habitat by increasing the number of cuts, and/or an effort to reduce 
the visual impact of clearcutting. 

The increase in the number of young, managed stand patches and the total 
amount of edge has implications for habitat potential of the landscapes. 
Although conclusive information on the effect of edges on vertebrates in 
western coniferous forest landscapes is not yet available, an increase in edge 
can benefit some species but prove detrimental to others (Yahner, 1988). For 
example, it has been observed that big game animals show an affinity for 
edges (Brown, 1985) and that some bird species occur more frequently on 
edges than in forest interiors (Rosenberg & Raphael, 1986). The fact that 
edge density increased in the present study areas over the past 15 years 
suggests that the habitat potential has increased for such species as elk, 
which can successfully utilize the edge environment. The increased dispersal 
of small clearcuts into the matrix of forest cover provides a corresponding 
increase in the amount of hiding cover close to forage areas used by the elk 
(Brown, 1985). 

Conversely, a number of other vertebrate species, such as the northern 
spotted owl, Townsend’s warbler Dendroica townsendi, and pileated 
woodpecker Dryocopus pileatus, may avoid edges (Bull, 1975; Brown, 1985; 
Rosenberg & Raphael, 1986). The increase in edge density indicates that 
habitat conditions for such species favoring interior forest have probably 
declined markedly. Specifically, the decrease in the mean distance of matrix 
to managed patch (as measured by GISfrag) and interior forest area are 
evidence of the decline in interior species habitat conditions. The manner 
and degree to which the decline affects interior species populations in the 
study areas are difficult to assess. Although Rosenberg and Raphael (1986) 
did not find a strong response of vertebrate communities in northern 
Californian landscapes which sustained a mean percent clearcut of 18% 
(roughly the same as the present study), the authors warned that the 
fragmentation in the region is a relatively recent phenomenon. Long-term 
vertebrate responses were not yet discernible. 

The continued use of dispersed clearcutting increases fragmentation of 
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forest landscapes at a rate more rapid than the rate of cutting on a per area 
basis. Given that the present cutting patterns are decreasing potential 
interior habitat at a rapid rate, alternative cutting patterns should be 
considered to reduce the loss of interior habitat and retain the area of large 
forested patches. Alternative models which aggregate cutting (Franklin & 
Forman, 1987) are available and may not require altering current standards 
and guidelines. 


CONCLUSIONS 

In an attempt to describe forest fragmentation, five groups of statistics were 
employed: patch abundance, patch shape, patch size, patch spacing, and 
matrix characteristics. By comparing two sets of data representing two dates 
over a 15-year period, we found that patch abundance, patch spacing 
measures, and matrix characteristics were most useful in capturing the 
amount of forest fragmentation over time. Patch size and shape statistics 
contribute information on specific characteristics of the individual patches 
and may be useful for applications designed to study specific interior and 
edge habitats or for the prescription of new clearcuts. In addition, a GIS 
fragmentation index was developed which proved sensitive to both the 
quantity and spacing of patches. The GIS provides an automated method of 
quantifying forest fragmentation to aid in forest and wildlife management 
decisions, and a means by which field and image data may be used in concert. 

Current concerns over forest fragmentation are typically related to a 
landscape condition in which forest islands occur in a matrix of managed 
forest plantations. This study suggests that on many Forest Service lands in 
the Cascade Range this condition is not yet realized, in contrast to many 
privately owned landscapes in the Cascades in which the matrix is the 
harvested area and the patches are the natural forest. Consequently, in 
characterizing fragmentation in some landscapes, characteristics of the 
matrix (the unmanaged forest) may be of more interest than characteristics 
of the patches (the managed plantations). Where the matrix is of interest, the 
interior area, the total edge, and the mean distance to the nearest managed 
patch (GISfrag) will be the useful descriptors of fragmentation. Additional 
research is needed to document and substantiate the relationship between 
forest landscape pattern and the subsequent wildlife/ecosystem response. 
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